Time-reversal symmetry breaking surface states of d-wave superconductors induced 
by an additional order parameter with negative T c 
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Surface states of d a! 2_ !/ 2-wave superconductors are studied using the Ginzburg-Landau (GL) the- 
ory. For a [110] surface it has been known that the time-reversal symmetry (T) breaking surface 
state, (ei±is)-wave state, can occur if the bare transition temperature of the s-wave order parameter 
(OP) is positive. We show that even if this bare T c is negative, it is possible to break T because the 
coupling to the spontaneously generated magnetic field may induce the s-wave OP. The T-breaking 
state is favored when the GL parameter k is small. 



PACS numbers: 74.20.De, 74.20.Rp 

I. INTRODUCTION 

Superconducting (SC) states of high-T c cuprates are 
known to have d x 2_ y 2-w&ve symmetry*^ Since the pair 
wave function of such an unconventional SC state has 
strong angular dependence, the effects of the pres- 
ence of surfaces, impurities are different from those 
in conventional s-wave superconductors. For exam- 
ple, it is possible to break the time-reversal symmetry 
(T) near a surface or a Joscphson junction by induc- 
ing the second component of the SC order parameter 
(OP)2£4£&LS&i&lLl&i&14 with a nontrivial phase dif- 
ference between the two OPs. In the case of a Josephson 
junction it may occur when the surface has [110] orien- 
tation, because the second SCOP induced by the tun- 
neling process can have phase difference ±7r/2 leading to 
a T-breaking state^S For a [110] surface faced to a vac- 
uum the necessary condition to break T seems to be that 
the bare transition temperature (T c ) of the second OP is 
positive^i ^ 2 .^ 

In this paper we examine the possibility to have a T- 
breaking surface state near the [110] surface of a d x 2_ y 2- 
wave superconductor when the bare T c of the additional 
OP is negative, namely, the second OP will not occur 
in the bulk even at zero temperature. We take an s- 
wave SCOP as the second component, since d x 2 _ y 2-w&ve 
and extended s-wave symmetries are natural candidates 
for superconducting states in the models with nearest- 
neighbor interactions (e.g., the t — J model). We will 
show that this kind of T-violation is possible, and that 
both the SCOPs and the magnetic field (vector poten- 
tial) should be treated self-consistently in order to de- 
scribe this situation correctly. It also turns out that the 
T violation may occur for a relatively small GL param- 
eter k (i.e., of the order of 10), when T c of the second 
OP is negative. Then the present mechanism may not 
be relevant to the T-violation in hole-doped cuprates in 
which k ~ 100. However, we expect the surface states of 
electron-doped cuprates may be described by the present 
theory, because some of the latter systems have much 
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II. GINZBURG-LANDAU EQUATION 

We consider a superconductor with tetragonal symme- 
try and assume only a <i :E 2_ y 2-wave SCOP, A d , is present 
in the bulk. An s-wave SCOP, A s , is taken into account 
as a possible second component when A^ is suppressed 
near the surface. For such a system the Ginzburg-Landau 
(GL) free energy is given as^ 



T = / dr [E (^|A M | 2 + ^|A M | 4 + iyDA,| 2 ) 

+ 7 i|aA£| 2 +72{A 2 (A:) 2 + (A3) 2 A 2 } 

+ K ds {(D x A d )(D x A s )* - (D y A d ){D y A s )* + c.c} 



where A is the vector potential and D = V — (27ri/$o)A 
is the gauge invariant gradient with $o = hc/2e being 
the magnetic flux quantum. Coefficients a M (oc T — T CM ), 
/3 M , Kfj,, 71, 72 and K ds are real, and we assume T cd > 0, 
while T cs can be both positive and negative. The 72 
is one of the terms which determine the relative phase 
of OPs, (j) ds {= 4> d - (p s ;A^ = |A M |exp(i</> M )). We take 
72 > 0, because this choice would lead to the (d ± is)- 
state (<p ds = ±7r/2) instead of the (d ± s)-state (<f) ds — 
0, 71"). In the former case the nodes of the (i-wave state 
are removed and the more condensation energy can be 
gained. It is also to be noted that 71 ± 272 is positive in 
usual weak-coupling model, since two OPs compete each 
other. Now we rewrite T in the dimensionless unrfeii to 
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see the parameter dependence of the model more clearly, 
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where 77^ = A^/Aq (/U = d, s) with Aq 
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being the bulk ci-wave OP. r was rescaled using the co- 
herence length for the d-wave OP, ^ (= y/IQj\aj\), as 
r — > r/^d, and D = V — ia/rt. Here a = A/(*/2H c t;d), 
and the magneti c field is measured in units of ^/2H C , 
where H c = \J AnrajTfid is the thermodynamic criti- 
cal field, k = Xd/S,d is the GL parameter with = 
cj)o/{2y/2TTH c ^d) being the penetration depth for the bulk 
d-wave superconductor. The parameters in Eq.(2) are 
defined as a s = a s /\a d \, f3 s = f3 s /f3 d , K s = K s /K d , 
7i = 71/A2, 72 = 12/Pd and K ds = K ds /K d . 

Usually the surface effect is described by the 
second-order surface GL free energy, J- S i = 

dS OnuVfliJIui where integration is carried 

sf fi,v=d,s 

out on the surface. Using the symmetry argument we 




FIG. 1: Schematic of a [110] surface of a d a .2_ !/ 2-wave super- 
conductor with tetragonal symmetry, x and y are parallel to 
the crystal a and b axes, respectively. 



find gds — 9sd — 9o cos 29 where 9 is the angle between 
the surface and the crystal a-axis with go being a 
constant. This term could also determine <pds, and it 
leads to the (d ± s)-state in the case of a [100] surface 
(9 = 0), since the 72 term is higher order than the gds 
term. However, gd s vanishes for a [110] surface (9 = 45°) 
which we consider in the following. The g^ term will 
represent the suppression of 77^ near the surface. Instead 
of using gdd we impose the condition rjd — at the [110] 
surface, because the d x 2_ y 2-w&ve SCOP should vanish 
there. Since the s-wave SCOP is only little affected 
by the presence of the surface, we take g ss — 0. (In 
numerical calculations we have checked that taking small 
positive g ss will not change the results qualitatively.) 



In order to consider the [110] surface we transform the 
coordinate system, (x,y,z) — * (x, y, z). Here x (y) is 
parallel to the crystal a (b) axis (z is parallel to the 
surface) , and i and y axes arc perpendicular and parallel 
to the surface, respectively. (See Fig.l.) In the free 
energy density only the K ds term is changed under this 
transformation to 



2Kd 



-aylm(r]*d x r)d ~ Vdd x r]* s ) 



(3) 



where we have assumed that the system is uniform along 
the surface, and the gauge freedom was taken as a = 

CLy(x)ey. 

The expression for the supercurrent is obtained by 
varying the electronic part of T (i.e., except the last 
term) with respect to a. Since the surface is faced to the 
vacuum, the x component, J x , should obviously vanish. 
(We have numerically checked that J x actually vanishes.) 
The y component, Jy, and that in the the dimensionlcss 
unit, are given as 



jy " Jy/ ( 
1 

4?r 
K d 



V2H C 
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ay(\r, d \ 2 + K s \rj s \ 2 ) 



(4) 



III. 



lm(rj*d x r]d - Vdd x r]*) 

K 



SURFACE STATE AND SPONTANEOUS 
CURRENT 



We numerically solve the problem by employing the 
quasi-Newton method^ to minimize the free energy T 
under the condition rj d {x = 0) = 0. We minimize T 
with respect to all variables, i.e., rjd, r] s and ay. Note 
that the Maxwell's equation is taken into account in this 
procedure, and we call this as "fully self-consistent cal- 
culation" . For the sake of comparison we will also show 
the results by treating only rjd and rj s self-consistently. 

First let us consider the case of a s < (i.e., T < T cs ). 
In this case, we would get finite rj s if rjd were absent. 
However, for T c d > T cs the stability condition for rj s in 
the bulk is given as, a s + (71 — 272) Iri^ 2 < 0, so the tran- 
sition temperature of rj s is lower than the bare one, T cs , 
and 7y s would be totally suppressed if T c d 3> T cs . Near 
the surface or impurities the situation can be different. 
There 77^ may be finite because the dominant SCOP, 77^, is 
suppressed. In Fig. 2 the spatial variations of the SCOPs 
near the surface are shown. ?7 S gets finite near the sur- 
face while rid is suppressed. The relative phase (j>ds will 
be determined by 72 and Kd s terms, and the former fa- 
vors cj>ds = ±7r/2 as mentioned. From Eq.(3) we see that 
the Kds term also favors 4>ds = ±tt/2, and ay will be 
spontaneously generated. (We take 77^ to be real and 
ay = in the bulk, i.e., x — + 00.) Numerical calculations 
show that rjd is real for all x, and that (f>ds = ±7r/2 where 
rjs is finite. This indicates that a T-violating (d + is)- 
wave surface state with a spontaneous magnetic field b z 
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(= dxdy) and a supercurrent jy occurs near the surface. 
The spatial distributions of b z and jy are presented in 
Fig.3. 

In order to see the role played by the vector potential, 
we investigate the same problem by setting a.y — ev- 
erywhere. Namely we treat only SCOPs self-consistcntly. 
When a,y is set to zero, the spontaneous current jy has 
contributions from only the spatial variations of SCOPs 
(i.e., the last line of Eq.(4)), and we calculate the mag- 
netic field from jy using Maxwell's equation, jy(x) — 
_ i db^X For - Q th regults for the SCOPs look 

similar as in the fully self-consistent calculations. The T- 
breaking (d + is)-state occurs as shown in Fig. 2. On the 
contrary, the behaviors of b z and jy are different in that 
jy always has the same sign, and that b z is a monotonous 
function of x. These results are not correct even qual- 
itatively as well as in a quantitative sense. Integration 
of the Maxwell's equation with the boundary condition 
6 z (±oo) = leads to /_ dxjy(x) = 0, implying that the 
averaged current should vanish.- This is the case for the 
fully self-consistent calculation but not in the case where 
the magnetic field is not treated self-consistently, because 
of the absence of the screening effect in the latter. 
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FIG. 3: Spatial variations of b z and jy. Parameters used 
are the same as in Fig. 2. (a) j y and (b) b z in the fully self- 
consistent calculation, (c) jy and (d) b z in the simplified one 
without treating ay self-consistently. Note b z and j y are in 
the dimensionless unit. 



energy gain is solely coming from the iCj s term. The 
spatial variations of b z and jy are shown in Fig. 5. 
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FIG. 2: Spatial variations of SCOPs for a s = -0.2, j3 s = 
0.2, K s = 0.5, 71 = 0.5, 72 = 0.1, K ds = 0.3 and k = 16. Note 
that all SCOPs are normalized by the bulk d-wave OP, and 
x — corresponds to the surface faced to the vacuum, (a) 
Ker/d and (b) Im?7 s in the fully self-consistent calculation, (c) 
Rer/d and (d) ImT) s in the simplified one without treating ay 
self-consistently. 

Next we consider the case of a s > 0, i.e., T > T cs . 
Note that T cs may be negative, in which case rj s will not 
occur in the bulk at T = even when r\d is absent. The 
results for the SCOPs are depicted in Fig. 4. (Here the 
GL parameter is taken to be k = 16.) It is seen that fi- 
nite Im(?7 s ) is obtained, though we naively expect rj s = 0. 
This is because the iCj s term couples SjRe^) bilinearly 
to aylm(ri s ). It may induce the state with Im(ry s ) ^ and 
b z 7^ 0, but the state with rj s = and b z = may also be 
a self-consistent solution. Numerical calculations show 
that the former one has the lower energy, and thus the 
time-reversal symmetry is violated spontaneously. Here 
|cfc a |, (3 S and K s were taken to be much smaller than those 
in Fig. 2. Otherwise the T-violation will not occur, be- 
cause these terms cost the energy for a s > and the 



FIG. 4: Spatial variations of SCOPs for a 3 = 0.01, f3 s = 
0.01, K s = 0.04, 71 = 0.5,72 = 0.1, K ds = 1.0 and k = 16. 
(a) Kerjd and (b) Imr7 s in the fully self-consistent calcula- 
tion, (c) Kerjd in the simplified one without treating ay self- 
consistently. 



In the case of a s > 0, the results with or without treat- 
ing the vector potential self-consistently are completely 
different. If we do not take into account the a,y term, 
r] s will never appear, since there is no mechanism to de- 
rive finite r/ s . Thus neither the spontaneous current nor 
the spontaneous field can occur. It implies that the T- 
violation near the surface cannot be described in this kind 
of simplified treatment for the superconductors in which 
the second SCOP has negative T c . 

In order to see the dependence on k we show the results 
for a larger k (k= 19) in Fig. 6 and 7. It is seen that \r/ s \, 
\b z \ and \jy\ are much smaller than those for k = 16. 
This k dependence can be understood as follows. r\d is 
suppressed in the region near the surface (x < £d), and 
rj s and ay would be finite there if T is broken. On the 
other hand the magnetic field b z would be finite in the 
region x < A<j. When k is large, the loss of energy due 
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FIG. 5: Spatial variations of b z and jy. Parameters used 
are the same as in Fig. 4. (a) jy and (b) b z in the fully self- 
consistent calculation. 



to finite b z in the large region < x < A^) overwhelms 
the energy gain coming from the Kd s term which acts 
only in the small region x < £4. Thus for large k the 
T-violation is not favored. If the larger value of Kd s is 
taken, the T-breaking state can occur for larger k. But 
the natural assumption seems to be Kds < Kd (Kds < 1), 
so that the T-violation may occur for k of the order of 
10. (On the contrary the T-violation may occur for much 
larger k in the case of a s < 0, because the energy can 
be gained by not only Kd s but also a s term.) It implies 
that the present mechanism may not be relevant to hole- 
doped cuprates in which k ~ 100, but it may describe 
the surface states of electron-doped cuprates which have 
smaller k. 




FIG. 6: Spatial variations of SCOPs. Parameters used are 
the same as in Fig. 4 except k = 19. (a) Rer^ and (b) Imr; 3 in 
the fully self-consistent calculation, (c) Ker/d in the simplified 
one without treating dy self-consistently. 



If we assume H c =1T, the maximum values of \B Z | and 
|Jji| are 2.5 x 10 _1 T and 3.7 x 10A/cm 2 , respectively, for 



k = 16. For k — 19 they are 8.6 x 10~ 2 T and 1.2 x 
1 OA/cm 2 , respectively. These values rapidly decrease as 
k increases, and the T-breaking state disappears as n 
exceeds 19 for the parameters used here. If we compare 
these values with experiments, it should be noted that 
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FIG. 7: Spatial variations of b z and jy. Parameters used 
are the same as in Fig. 6. (a) jy and (b) b z in the fully self- 
consistent calculation. 



surface roughness will reduce \B Z \ and \J y \, because T- 
violation is most favored in the case of 9 = 45° i 2 - (When 
9 ^ 45°, gds will be finite and the T-violation is not 
favored.) 



IV. SUMMARY 

We have examined the role played by the vector po- 
tential concerning the occurrence of surface states with 
spontaneously broken T in d a ,2_j,2-wave superconductors. 
It has been known that the T-breaking state may natu- 
rally appear if the bare T c of the additional OP is pos- 
itive. For the Josephson junction composed of d x 2_ y 2- 
wave and other superconductors, tunneling may induce 
second component of SCOP and thus T may be broken. 
In these cases the T-breaking states may be described 
without treating the vector potential self-consistently. In 
this paper it was shown that the surface state of a d x i_ y i- 
wave superconductor may break T even when the bare 
T c of the second SCOP is negative. However, to describe 
this situation correctly not only the SCOPs but also the 
vector potential must be treated on an equal footing. 
In the present mechanism the T-violation may occur for 
rather small values of the GL parameter k (< 20), so that 
it may not be relevant to hole-doped cuprates. We expect 
the present theory may be used to describe the surface 
states of electron-doped high-T c cuprates, because their 
k are much smaller than those of hole-doped systems. 
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